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ABSTRACT 

Down-looking  ground-penetrating  radar  (DLGPR)  has  been  used  extensively  for  buried  target  detection.  For  operational 
implementations,  the  sensor  is  used  in  direct-detection  mode,  where  algorithms  process  data  while  the  system  moves 
down  roadways.  Decisions  are  made  before  a  system  passes  over  the  target.  Change  detection  works  by  passing  over  an 
area  before  and  after  targets  are  buried.  By  comparing  before-and-after  data,  change  detection  can  improve  DLGPR 
performance,  but  it  also  has  inherent  operational  limitations.  Performance  enhancements  include  mitigating  the  effects  of 
anomalies  not  associated  with  targets  and  increasing  the  detection  probabilities  of  deeper  targets  through  indirect  means. 
In  the  latter  case,  deeply  buried  targets  that  do  not  appear  in  the  GPR  data  can  be  indirectly  detected  using  change- 
detection  methods  if  the  patch  of  ground  where  the  target  is  buried  has  been  significantly  modified  from  its  original 
undisturbed  state.  In  this  paper,  we  explore  decision-based  change-detection  approaches  for  enhancing  the  performance 
of  a  DLGPR  system  and  enumerate  the  limitations  of  the  approach. 

Keywords:  change  detection,  down-looking  ground-penetrating  radar,  DLGPR,  mine  detection 

1.  INTRODUCTION 

Most  of  the  down-looking  ground  penetrating  radar  (DLGPR)  systems  employed  for  buried  target  detection  have 
concentrated  on  the  problem  of  directly  detecting  the  anomaly  in  the  GPR  data  caused  by  the  buried  object.  Furthermore, 
most  DLGPR  algorithms  attempt  to  reduce  false-alarm  rates  through  some  form  of  discrimination,  in  which  features  are 
extracted  and  classifiers  are  trained  on  known  targets  and  false  alarmsk  With  change  detection,  feature  extraction  and 
classification  are  de-emphasized,  the  reductions  in  false  alarms  resulting  instead  from  comparisons  of  the  ground  before 
and  after  targets  are  buried. 

For  DLGPR  systems,  the  largest  response  is  typically  due  to  the  radar  reflecting  off  the  ground.  Most  DLGPR  algorithms 
remove  the  ground-bounce  response  as  a  first  pre-processing  step.  The  remaining  subsurface  response  is  then  used  to 
detect  buried  targets.  DLGPR  systems  are  not  able  to  detect  some  targets  due  to  the  soil  conditions  in  which  they  are 
buried,  the  depth  at  which  they  are  buried,  or  both.  Those  targets  simply  do  not  show  up  in  the  subsurface  response  of  the 
radar.  But  when  the  ground  is  scanned  before  the  targets  are  buried,  change  detection  can  be  used  to  indirectly  detect 
targets  that  otherwise  cannot  be  detected  directly.  Despite  some  operational  limitations  imposed  on  the  implementation 
of  change  detection,  the  approach  can  still  be  assessed  for  its  efficacy.  In  this  paper,  we  apply  various  algorithms  to  data 
collected  in  support  of  a  change-detection  test.  Sample  GPR  data  associated  with  algorithm  outputs  sheds  light  on  the 
process.  We  provide  receiver  operator  characteristic  (ROC)  curves  that  demonstrate  enhanced  performance  using  change 
detection. 


2.  LIMITATIONS  OF  CHANGE  DETECTION 

The  change-detection  results  and  conclusions  presented  in  this  paper  are  encouraging  but  have  many  limitations  that 
must  be  overcome  before  these  techniques  can  be  used  in  the  field.  Limitations  of  implementing  change  detection  for 
DLGPR  systems  include  a  reliance  on  differential  Global  Positioning  System  (DGPS)  to  accurately  correlate  alarms 
from  one  excursion  to  another,  the  misidentification  of  changes  in  the  soil,  and  logistical  implementation  hurdles. 

The  results  presented  in  this  paper  are  derived  from  data  collected  using  real-time  kinematic  (RTK)  DGPS,  which  is 
accurate  to  a  few  centimeters.  RTK-DGPS  requires  the  use  of  a  base  station  at  known,  surveyed  locations,  with  radios 
transmitting  differential  corrections  from  the  base  station  to  a  receiver  on  the  moving  DLGPR  system.  The  resulting 
accuracy  of  the  alarm  positions  allows  for  the  correct  association  of  alarms  caused  by  the  same  subsurface  source — from 
one  excursion  to  another  or  from  one  day  to  another.  The  required  accuracy  of  the  alarm  positions  is  in  part  determined 
by  the  size  of  the  targets  that  the  system  is  designed  to  detect.  Since  the  targets  are  typically  less  than  1/2  meter  in  size, 
linking  distances  are  approximately  the  same  length.  If  a  GPS  system  without  differential  corrections  was  used  to 


compute  alarm  locations  with  an  accuracy  of  a  meter  or  more,  it  would  be  difficult  to  accurately  link  an  alarm  due  to  a 
subsurface  object  from  one  excursion  to  another.  At  such  accuracy,  an  alarm  from  a  subsurface  object  could  be  reported 
as  far  away  as  several  meters  from  an  alarm  from  the  same  subsurface  object  generated  from  a  prior  excursion.  Although 
some  types  of  navigation  systems  can  enhance  the  accuracy  of  GPS  systems  that  do  not  benefit  from  differential 
corrections,  these  are  not  considered  here. 

Alarms  produced  by  automatic  target  recognition  (ATR)  algorithms  on  board  the  DLGPR  system  depend  on  soil 
conditions.  It  is  therefore  expected  that  change  detection  using  DLGPR  systems  will  suffer  from  high  false-alarm  rates 
when  soil  conditions  vary  from  one  excursion  to  another.  For  example,  if  the  soil  is  dry  one  day  and  moist  the  next,  the 
alarms  generated  by  the  ATR  can  differ  significantly  in  number  and  in  confidence  value.  As  a  result,  many  alarms  will 
not  have  a  natural  match  from  one  excursion  to  another.  Change  detection  requires  at  least  a  match  in  alarm  locations 
and  sometimes  requires  that  alarm  confidence  values  be  nearly  the  same.  In  conditions  in  which  roadways  are  traveled 
by  vehicular  or  foot  traffic,  it  is  expected  that  the  impact  on  the  soil  conditions  of  the  road  surface  and  shallow 
subsurface  will  change  from  day  to  day,  or  from  hour  to  hour.  These  changes  can  potentially  give  rise  to  alarms  even 
though  they  have  nothing  to  do  with  buried  targets. 

To  implement  change  detection  into  roadway-monitoring  scenarios,  the  alarms  from  all  excursions  must  be  accurately 
shared  between  each  of  the  DLGPR  vehicles.  This  requirement  imposes  a  strict  set  of  protocols  to  be  followed.  The 
effectiveness  of  change  detection  relies  on  the  continual  monitoring  of  a  given  roadway,  which  places  additional 
logistical  burdens  on  those  responsible  for  monitoring  the  road. 


3.  DLGPR  DATA  FUNDAMENTALS 

To  evaluate  various  approaches  to  direct  detection  and  change  detection,  we  study  the  sensor  data  in  their  fundamental 
forms,  adopting  the  terminology  used  in  Daniels  book  on  GPR^.  The  fundamental  response  of  any  DLGPR  system  is  the 
A-scan,  which  is  associated  with  a  particular  down-track  location  and  a  particular  across-track  channel.  When  we 
examine  DLGPR  data  from  most  systems,  we  typically  differentiate  raw  A-scans  from  processed  A-scans.  Figure  1 
shows  an  A-scan  response  collected  over  a  buried  metallic  target.  The  largest  response  is  from  the  ground,  followed  by 
two  distinct  high-magnitude  reflections  from  the  target  itself.  Nearly  all  other  parts  of  the  response  are  zero,  including 
the  response  before  the  ground  bounce,  the  response  between  the  ground  and  the  buried  target,  and  the  response  after  (or 
deeper)  than  the  buried  target. 
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Figure  1.  DLGPR  A-scan  for  a  buried  metallic  target. 

If  we  rotate  the  A-scan  of  figure  I  90  degrees  so  that  the  ground  bounce  is  at  the  top  and  then  plot  neighboring  A-scans 
to  the  left  and  right  of  it,  we  get  the  resultant  waterfall  plot  shown  in  figure  2.  There  are  51  across-track  channels  in  this 
particular  DLGPR  system.  All  the  A-scans  of  figure  2  come  from  channel  40,  which  happens  to  be  the  channel  in  which 
the  target  is  centered.  Thus,  the  adjacent  A-scans  are  from  successive  along-track  scans,  which  are  spaced  approximately 


every  0.05  meters.  In  figure  2,  the  x-axis  is  the  local  scan  number,  and  so  the  51  A-scans  shown  were  collected  over 
about  2.5  meters  of  along-track  distance.  (We  use  the  terms  along-track  and  down-track  interchangeably). 


Figure  2.  DLGPR  waterfall  plot  of  A-seans  for  a  buried  metallie  target. 

If  instead  of  plotting  x-y  A-scans,  we  create  a  grayscale  image  of  the  GPR  responses,  we  get  an  along-track  B-scan,  as 
shown  in  figure  3.  Here,  the  peaks  and  valleys  of  the  GPR  responses  have  been  mapped  to  a  color  on  the  gray-scale  color 
bar,  with  white  being  a  high-magnitude  peak  and  black  being  a  high-magnitude  valley.  We  observe  that  the  metallic 
target  manifests  itself  as  a  parabola  in  this  B-scan  view  of  the  GPR  data. 
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Figure  3.  DLGPR  along-track  B-scan  for  a  buried  metallic  target. 

We  now  construct  a  two-dimensional  C-scan  from  the  three-dimensional  volumetric  GPR  data.  The  matrix  dimensions  of 
the  C-scan  are  given  by  the  number  of  down-track  scans  and  by  the  number  of  across-track  channels.  In  this  way,  the  C- 
scan  can  be  thought  of  as  a  bird’s-eye  view  of  the  lane,  or  roadway,  which  was  scanned  by  the  DLGPR  system.  The 
simplest  way  to  generate  a  C-scan  is  to  take  a  slice  of  the  three-dimensional  GPR  data  at  a  fixed  depth,  or  time  bin. 
Alternatively,  there  are  countless  ways  to  transform  the  three-dimensional  GPR  data  to  a  two-dimensional  C-scan^.  A 
common  technique  is  to  take  the  sum  of  the  squares  of  each  A-scan.  Other  methods  include  KSum"^,  Disturbed  Soil 
Metric^,  and  LMS^.  Figure  4  shows  a  C-scan  corresponding  to  the  same  metallic  target  shown  in  the  B-scan  of  figure  3, 


in  this  case  derived  from  a  variant  of  the  LMS  algorithm.  Essentially,  the  volumetric  three-dimensional  GPR  data  are 
processed  and  collapse  to  two  dimensions.  Targets  manifest  themselves  as  blobs  in  the  C-scan  images,  where  white 
pixels  are  of  the  highest  magnitude  and  black  pixels  are  of  the  lowest  magnitude. 
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Figure  4.  DLGPR  C-scan  for  a  huried  metallic  target. 

An  illustrative  way  of  presenting  the  GPR  data  at  a  target  or  alarm  location  is  to  show  the  down-track  B-scan,  across- 
track  B-scan,  and  a  C-scan  projected  onto  the  faces  of  a  cube.  Figure  5  shows  a  DLGPR  cube  at  the  location  of  the 
buried  metallic  target.  High-intensity  blobs  in  the  C-scan  are  mapped  to  a  color  scale  where  magenta  represents  the 
highest  magnitude  response. 
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Figure  5.  DLGPR  cube  for  a  buried  metallic  target. 


4.  DECISION-BASED  CHANGE  DETECTION 


4.1  General  Approach 

In  decision-based  change  detection,  alarms  are  compared  from  excursions  before  and  after  targets  are  buried.  A 
comparison  of  their  locations  is  all  that  is  necessary  to  gain  a  benefit  from  this  form  of  change  detection.  But  if  we  also 
use  the  confidence  levels  of  the  alarms,  we  can  improve  the  change-detection  result. 

Figure  6  depicts  the  operational  and  testing  scenario  that  suggests  a  change-detection  approach  for  detecting  buried 
targets.  Each  of  the  sub-figures  represents  an  event  in  time.  In  figure  6a,  a  DLGPR  system  scans  a  roadway  before  any 
targets  are  buried.  The  system  generates  a  list  of  suspected  target  locations,  or  alarms,  which  are  shown  as  numbers.  The 
essential  elements  of  an  alarm  are  a  unique  alarm  number,  a  GPS  location,  and  an  ATR  scalar  output — also  referred  to  as 
a  confidence  level.  The  higher  the  ATR  output,  the  greater  the  confidence  that  the  alarm  is  due  to  a  target  and  not  a  false 
alarm.  Since  no  targets  were  present  in  the  roadway  when  the  system  made  its  first  pass,  the  four  alarms  in  figure  6a  are, 
by  default,  false  alarms.  For  this  analysis,  we  are  not  interested  in  the  source  of  these  false  alarms.  Since  no  known 
targets  were  encountered,  the  probability  of  detection  (Pd)  is  not  defined. 
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Figure  6.  Depiction  of  change  detection  approach. 

Between  the  scenarios  depicted  in  figures  6a  and  6b,  four  targets  are  buried  in  the  roadway.  Targets  are  depicted  as  gray 
circles,  and  the  halo  radii  are  depicted  as  dashed  gray  circles  surrounding  the  targets.  Halo  radii  are  the  allowable  miss 
distances  between  an  alarm  and  a  target.  Alarms  within  the  halo  of  a  target  are  deemed  detections,  and  alarms  outside  the 
halos  are  deemed  false  alarms.  In  figure  6b,  we  see  the  result  of  a  second  pass  by  the  DLGPR  system.  A  new  set  of 
alarms  is  generated,  this  time  represented  by  numbers  in  red.  After  targets  are  buried,  a  total  of  six  alarms  are  generated 
by  the  direct  detection  algorithm  of  the  DLGPR  system.  We  see  in  figure  6b  that  alarms  1,  4,  and  5  fall  within  the  halo  of 
three  of  the  four  targets,  but  alarms  2,  3,  and  6  fall  outside  the  halos  of  the  targets.  Target  2  is  not  detected  (it  could  be  a 
deeply  buried  target  that  the  system  was  not  designed  to  detect).  The  direct-detection  result  is  three  out  of  four  targets 
detected,  with  three  false  alarms. 


4.2  Change  Detection  Using  Only  Locations  of  Alarms 

We  now  apply  change  detection  by  overlaying  the  alarms  from  the  excursion  before  targets  were  buried  (black  numbers) 
with  the  alarms  from  the  excursion  after  targets  were  buried  (red  numbers).  The  result  is  shown  in  figure  6c.  For  our  first 
attempt  at  change  detection,  we  compare  alarm  locations  from  the  second  excursion  with  alarm  locations  from  the  first 
excursion  and  remove  alarms  that  are  spatially  correlated.  For  our  analysis,  we  used  a  linking  distance  of  0.5  m.  In  figure 
6c,  red  circles  are  placed  around  alarm  pairs  that  meet  the  linking-distance  criterion.  After  eliminating  spatially 
correlated  alarms,  we  are  left  with  fewer  alarms  for  the  second  excursion,  as  shown  by  the  red  numbers  in  figure  6d.  For 
our  initial  approach,  we  do  not  retain  any  of  the  alarms  from  the  first  excursion  (black  numbers).  With  reference  to  figure 
6d,  we  see  that  there  are  only  three  remaining  alarms.  The  change-detection  result  is  two  out  of  four  targets  detected, 
with  one  false  alarm.  For  this  case,  we  eliminated  two  of  the  three  false  alarms,  but  also  removed  one  of  the  three 
detections  of  targets.  It  is  the  goal  of  change  detection  to  reduce  the  false-alarm  rate  significantly,  with  only  modest 
reductions  in  the  probability  of  detection.  Optimally,  we  would  like  to  eliminate  all  the  false  alarms  without  sacrificing 
any  detections. 

4.3  Enhancing  Change  Detection  Performance  by  Using  Alarm  Confidence  Values 

We  reexamine  what  happened  with  target  3.  We  see  that  it  was  detected  in  figure  6b,  but  the  alarm  from  the  second 
excursion  (red  alarm  number  4)  was  then  eliminated  by  virtue  of  being  near  an  alarm  from  the  first  excursion  (black 
alarm  number  4).  We  might  say  that  alarm  number  4  from  the  first  excursion  (black  number)  was  in  the  wrong  place  at 
the  wrong  time.  How  might  we  modify  our  change-detection  algorithm  so  that  red  alarm  number  4  is  not  eliminated? 
What  additional  information  is  at  our  disposal,  other  than  the  location  of  the  alarms?  The  answer  is  that  we  can  compare 
not  only  the  location  of  first  and  second  excursion  alarms,  but  also  their  direct-detection  algorithm  confidence  values.  If 
the  confidence  values  of  spatially  correlated  alarms  are  similar,  we  proceed  to  eliminate  the  second  excursion  alarms.  By 
doing  this,  we  are  essentially  saying  that  the  ground  has  not  changed  from  the  first  excursion  to  the  second  excursion. 
Whatever  subsurface  anomaly  gave  rise  to  a  direct-detection  alarm,  it  appears  essentially  the  same  from  the  first  pass  to 
the  second  pass.  That  is,  these  alarms  are  deemed  as  arising  from  the  same  source,  which  means  that  no  change  occurred 
in  the  roadway  at  these  locations.  In  the  case  that  the  confidence  value  of  the  second  excursion  alarm  is  greater  than  the 
confidence  value  of  the  first  excursion  by  a  certain  percentage,  we  do  not  remove  the  second-pass  alarm.  Here,  the 
assumption  is  that  something  changed  in  the  ground  to  give  rise  to  an  alarm  with  a  significantly  greater  confidence  value. 
For  our  problem  it  was  the  introduction,  or  burial,  of  a  target.  This  case  is  depicted  in  figure  7,  in  which  the  second 
excursion  alarm  (red  alarm  4)  has  a  confidence  value  of  9.2,  which  is  deemed  significantly  greater  than  the  first 
excursion  alarm  (black  alarm  4),  which  has  a  confidence  value  of  2.1.  Figure  8  shows  the  results  of  using  this  enhanced 
decision-based  change-detection  algorithm.  Now  the  number  of  false  alarms  has  been  reduced  to  one,  as  before,  and  the 
number  of  detections  is  restored  to  3. 


Figure  7.  Alarm  numbers  and  eonfidenees  assoeiated  with  target  4. 
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Figure  8.  Depietion  of  ehange  deteetion  result  when  alarm  eonfidenee  is  used,  in  addition  to  alarm  loeation. 


We  now  turn  to  the  case  of  target  2,  which  we  suppose  is  a  deeply  buried  target  that  the  DLGPR  cannot  detect  directly. 
Using  change  detection,  are  there  means  by  which  deeply  buried  targets  can  be  detected  indirectly  by  virtue  of 
comparisons  of  the  same  patch  of  ground  before  and  after  the  target  is  buried?  Does  the  GPR  data  hold  more  clues  or 
evidence  of  a  deeply  buried  target,  other  than  the  typical  direct-detection  approach  of  detecting  and  discriminating 
hyperbolic  anomalies?  In  the  next  section,  we  examine  one  such  method. 


5.  CHANGE  DETECTION  USING  GROUND-BOUNCE  METRIC 

A  metric  was  developed  to  identify  a  specific  phenomenon  that  sometimes  occurs  when  targets  are  buried.  The 
phenomenon  can  be  particularly  pronounced  for  deeply  buried  targets  since  they  tend  to  require  a  larger  hole  to  be  dug 
than  smaller  targets.  The  metric  allows  the  system  to  indirectly  detect  deeply  buried  targets,  which  likely  cannot  be 
directly  detected  by  the  GPR.  The  metric  can  be  implemented  as  a  change-detection  method,  where  the  current  view  of 
the  ground  is  examined  and  compared  to  the  ground  at  the  same  location  at  a  previous  time,  or  as  a  direct-detection 
method,  where  only  the  current  view  of  the  ground  is  analyzed. 

Figure  9  depicts  the  phenomenon  that  the  metric  was  created  to  detect.  The  set  of  B-scans  on  the  left  correspond  to  a 
location  traversed  by  a  DLGPR  without  a  target  present.  The  B-scans  on  the  right  correspond  to  the  same  location  after  a 
target  was  buried  at  a  deep  depth.  Physics  limitations  prevent  the  DLGPR  from  receiving  a  response  from  a  target  this 
deep,  and  thus  the  target  cannot  be  seen  in  the  B-scan  images.  But  there  is  evidence  of  the  target  burial  when  the  ground- 
bounce  location  is  examined.  On  the  left,  the  ground  is  relatively  flat  in  the  down-track  B-scan  and  almost  completely 
flat  in  the  across-track  B-scan.  After  the  target  is  buried,  a  depression  is  centered  at  the  target  location,  as  shown  in  the 
B-scans  on  the  right. 
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Figure  9.  Across-track  and  down-track  B-scans  of  a  location  before  and  after  target  burial. 

The  aptly  named  ground-bounce  metric  utilizes  ground-bounce  location  information  to  attempt  to  identify  the  scenario 
where  there  is  a  depression  at  the  location  of  a  deeply  buried  target.  Unlike  the  previous  decision-based  change-detection 
methods,  this  metric  is  a  stand-alone  routine  applied  to  three-dimensional  GPR  data  to  create  alarm  files  that  can  be 
fused  with  other  pre-screener  alarm  files  in  hopes  of  attaining  a  better  overall  ROC  curve  due  to  the  additional  detections 
of  deeply  buried  targets. 


The  algorithm  begins  by  calculating  the  ground-bounce  location  for  every  channel  at  each  down-track  scan,  creating  a 
two-dimensional  C-scan  map  of  these  values.  Here,  the  ground-bounce  location  is  defined  as  the  time  bin  number 
corresponding  to  the  maximum  radar  amplitude  response  due  to  the  ground  bounce.  If  the  algorithm  is  being  utilized  as  a 
change-detection  method,  then  C-scan  maps  of  ground-bounce  locations  are  created  for  the  GPR  data  before  and  after 
targets  are  buried.  To  account  for  differing  antenna  heights  at  the  same  down-track  location  before  and  after  targets  are 
buried,  the  algorithm  calculates  the  mean  ground-bounce  value  for  all  channels  in  each  down-track  scan.  The  algorithm 
then  subtracts  each  ground-bounce  location  in  the  C-scan  map  from  the  mean  ground-bounce  location  corresponding  to 
that  particular  down-track  scan.  If  the  algorithm  is  being  applied  as  a  change-detection  method,  the  C-scan  map  from 
before  the  target  is  buried  is  used  to  calculate  the  mean  ground-bounce  location,  and  the  C-scan  map  corresponding  to 
after  the  target  is  buried  is  subtracted  from  the  mean.  This  method  could  also  be  utilized  as  direct  detection  of  disturbed 
earth  by  only  examining  the  C-scan  corresponding  to  after  the  target  is  buried.  In  this  case,  the  algorithm  calculates  the 
mean  ground-bounce  location  of  all  channels  in  a  down-track  scan  and  then  subtracts  the  ground-bounce  values  at  every 
channel  in  the  same  scan  from  this  location. 

A  Prewitt  filter  is  applied  to  the  resultant  C-scan  to  eliminate  striping  and  utilize  edge-detection  features  to  enhance 
responses  due  to  depressions.  Figure  10  shows  C-scans  for  the  same  location  pictured  in  Figure  9  where  a  target  was 
buried  at  a  deep  depth.  The  C-scan  map  on  the  left  was  created  by  subtracting  the  ground-bounce  location  at  each 
channel  and  down-track  scan  from  the  mean  ground-bounce  location  at  that  down-track  scan.  The  resulting  C-scan 
contains  vertical  striping,  which  is  removed  by  applying  a  Prewitt  filter.  The  effects  of  applying  the  Prewitt  filter  can  be 
seen  in  the  C-scan  map  on  the  right;  the  striping  is  gone  and  the  location  of  the  recently  buried  target  stands  out  as  a 
bright  spot  compared  against  the  background. 


Figure  10.  C-scan  maps. 

Alarm  files  are  created  from  the  filtered  C-scan  map  by  applying  a  threshold  to  the  image  and  looking  for  areas  where 
five  or  more  adjacent  pixels  have  values  above  threshold.  Figure  11  is  a  flowchart  that  outlines  the  entire  procedure  by 
which  the  ground-bounce  metric  is  applied  to  the  GPR  data  to  create  alarm  files. 


Figure  11.  Flowchart  depicting  how  alarm  files  are  created  using  the  disturbed  soil  metric. 


6.  RESULTS 

Figure  12  gives  ROC  curves  for  the  change-detection  approaches  and  algorithms  portrayed  in  figures  6  through  8.  The 
part  of  the  ROC  curve  shown  in  figure  12  is  deemed  to  be  of  operational  interest.  If  the  false-alarm  rate  (FAR)  is  any 
higher,  it  is  expected  that  the  system  would  not  be  embraced  as  a  tool  for  target  detection  because  the  system  would  need 
to  stop  too  frequently  at  locations  where  no  targets  were  buried.  Error  bars  on  Pd  are  represented  by  fainter  colored  lines. 
We  see  that  when  using  only  the  location  of  the  alarms,  the  change-detection  approach  enhances  detection  performance 
at  low  false-alarm  rates.  For  a  FAR  of  0.001  m“^,  the  Pd  for  direct  detection  is  about  0.55  (black  curve),  while  for  change 


detection  it  is  about  0.65  (blue  curve).  But  at  a  higher  FAR  of  0.006  m“^,  the  Pd  for  change  detection  (using  only  location 
of  alarms)  is  lower  than  the  Pd  for  direct  detection.  The  benefit  of  eliminating  false  alarms  with  change  detection  is  offset 
by  a  reduction  in  Pd  at  a  FAR  of  0.006  m“^.  But  when  we  also  utilize  the  confidence  values  of  the  alarms  we  improve  the 
change-detection  result.  By  accounting  for  some  cases  in  which  an  alarm  was  in  the  wrong  place  at  the  wrong  time  (see 
alarm  number  4  in  figure  6a),  the  Pd  at  0.001  m“^  increases  from  0.65  to  about  0.70,  and  at  the  higher  FAR  of  0.006  m“^, 
no  degradation  in  Pd  is  seen.  The  value  of  change  detection  is  evident  when  we  compare  the  FARs  at  a  fixed  Pd.  Using 
only  direct  detection,  we  see  that  the  ROC  curve  reaches  its  final  plateau  (for  the  false-alarm  regime  shown  in  figure  12), 
or  its  maximum  Pd  of  about  0.70,  at  a  false-alarm  rate  of  0.004  m“^.  But  by  using  change  detection  (using  location  and 
confidence  values),  a  Pd  of  0.70  is  realized  at  a  false-alarm  rate  of  only  0.001  m“^.  The  number  of  false  alarms  is  reduced 
by  a  factor  of  4.  This  has  implications  for  fielded  DLGPR  systems.  If  the  limitations  of  using  change  detection  in  a 
fielded  system  can  be  overcome,  there  is  a  potential  for  delivering  detection  performance  as  much  lower  false-alarm 
rates,  such  that  system  thresholds  can  be  changed  accordingly,  and  operational  tempos  can  be  dramatically  increased. 


Figure  12.  Direct-detection  results  (black  curve),  change-detection  results  using  location-only  (blue  curve),  and  change- 
detection  results  using  location  and  confidence  values  of  alarms  (red  curve),  for  all  test  lanes. 

We  are  still  left  with  a  ROC  curve  that  plateaus  at  a  Pd  of  0.70.  This  performance  is  primarily  due  to  the  limitations  of 
the  GPR  sensor,  which  in  this  case  was  not  designed  to  detect  deeply  buried  targets.  By  employing  the  ground-bounce 
metric  described  in  section  5,  we  see  potential  for  the  indirect  detection  of  deeply  buried  targets.  Figure  13  shows  ROC 
curves  for  deep  targets,  limited  to  one  test  lane.  Here,  the  ground-bounce  metric  applied  in  a  change-detection  sense 
greatly  enhances  detection  of  deeply  buried  targets.  Direct  detection  results  in  a  Pd  of  zero  (black  curve  is  aligned  with 
the  x-axis),  and  application  of  the  ground-bounce  metric  results  in  the  detection  of  half  the  deeply  buried  targets.  It 
remains  to  be  seen  if  the  technique  is  robust  when  applied  to  the  larger  data  set.  The  data  used  for  the  ROC  curve  in 
figure  13  did  not  include  some  of  the  rougher  lanes,  which  could  induce  false  alarms  due  to  greater  variance  in  the 
ground-bounce  locations  from  run  to  run. 


1 


Direct 

Ground  Bounce  Metric 
Using  Change  Detection 


Figure  13.  ROC  curves  for  deep  targets,  limited  to  one  test  lane. 


7.  CONCLUSIONS 

Although  employing  change  detection  in  locations  where  high-accuracy  GPS  is  denied  is  fraught  with  technical 
challenges,  the  results  in  this  paper  indicate  that  applying  change  detection  to  DLGPR  systems  has  the  potential  to 
enhance  detection  performance.  By  utilizing  techniques  that  indirectly  detect  deeply  buried  targets,  detection 
performance  can  potentially  be  enhanced  against  the  entire  target  set. 


8.  FUTURE  WORK 

The  next  step  is  to  fuse  the  ground-bounce  metric  alarms  with  other  pre-screener  alarms  to  enhance  overall  detection 
performance.  The  goal  is  to  increase  the  Pd  plateau  of  figure  12  without  incurring  a  false-alarm  penalty  that  ultimately 
degrades  detection  performance.  One  idea  is  to  leverage  a  phenomenon  in  which  deeply  buried  targets  give  rise  to  a  sort 
of  blurred  effect  in  the  three-dimensional  GPR  data.  In  these  cases,  subsurface  clutter  is  destroyed  when  the  target  is 
buried  deeply,  and  thus  the  amount  of  subsurface  energy  actually  decreases  when  the  target  is  buried.  This  would  lead  to 
local  anomaly  nulls,  rather  than  peaks. 
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